arXiv: 1504.08119vl [nucl-th] 30 Apr 2015 


Rho meson form factors in a confining Nambu—Jona-Lasinio model 


Manuel E. Carrillo-Serrano/ Wolfgang Bentz,^ Ian C. Cloet,^ and Anthony W. Thomas^ 

^ CSSM and ARC Centre of Excellence for Particle Physics at the Tera-scale, 
Department of Physics, University of Adelaide, Adelaide SA 5005, Australia 
^Department of Physics, School of Science, Tokai University, 

4-1-1 Kitakaname, Hiratsuka-shi, Kanagawa 259-1292, Japan 
^Physics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 

Elastic electromagnetic form factors for the p"*" meson are calculated in a Nambu-Jona-Lasinio 
model which incorporates quark confinement through the use of the proper-time regularization scheme. 
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the lattice QCD calculations, with the exception of the quadrupole moment and corresponding form 
factor, which may be related to a lack of spherical symmetry on the lattice. 
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I. INTRODUCTION 

The structure of hadrons presents a remarkable chal¬ 
lenge to the theory of strong interactions - quantum chro¬ 
modynamics (QCD) - and a critical feature of a hadron’s 
structure is its distribution of charge and magnetization, 
which is empirically related to its electromagnetic form 
factors [1]. The direct calculation of hadron form factors 
using QCD is currently only possible through lattice QCD, 
albeit limited to the low to moderate region. However, 
to gain insight into the relevant dynamical mechanisms 
behind the observed structure it is useful to work with 
models that approximate key features of QCD. An impor¬ 
tant focus for this comparison are the meson form factors. 
Because of their short lifetimes [2] they present a unique 
challenge experimentally - making both lattice QCD and 
model calculations critical. The pion form factor has been 
successfully measured over a wide range of four momen¬ 
tum transfer, while the vector meson form factors have 
not had the same amount of experimental exploration. 
However, the BABAR collaboration has measured the 
cross-section for the reaction e+ -I- e“ —>■ p’*' + p~ [3], 
which has been analyzed to garner information on the 
p-meson form factors [4]. 

The p form factors, or equivalently the polarization am¬ 
plitudes, have been calculated using a variety of methods, 
for example, phenomenological models [5, 6], constituent 
quark models in the light front framework [7-14], QCD 
sum rules [15-18] and the Dyson-Schwinger equations [19- 
22]. The first attempts to compute p form factors using 
lattice QCD were reported in Refs. [23, 24] in the quenched 
framework. The recent work of Owen et al. [25] and Shultz 
et al. [26] give two independent lattice QCD calculations 
based upon different approaches. These lattice results, 
and the previous work with quark models, provides a solid 
background for comparison with results computed within 
other models. 

In this work we extend the p-meson form factor calcu¬ 
lation of Ref. [27], where the focus was a comparison with 
the axialvector diquark form factors which formed a criti¬ 


cal part of a nucleon form factor calculation. Here we use 
the same conhning version of the Nambu-Jona-Lasinio 
(NJL) model [28, 29] to investigate the quark mass de¬ 
pendence of the p form factors, and perform a detailed 
comparison with the lattice QCD results of Refs. [25, 26]. 
Following Ref. [27] we include the dressing of the quark- 
photon vertex from the inhomogenerous Bethe-Salpeter 
equation and a pion cloud, which are critical for a good 
agreement with lattice results. Similar finding were made 
in Ref. [30], where the same framework was applied to 
the TT and K form factors. 

The outline of the paper is as follows: In Sec. II we 
briefly review the NJL model as applied to qq bound 
states and the calculation of the p electromagnetic form 
factors is discussed in Sec. HI. The results are compared 
to those from lattice QCD and various quark models in 
Sec. IV and Sec. V presents our conclusions. 

II. NAMBU-JONA-LASINIO MODEL 

In its original formulation the NJL model successfully 
encapsulated the effects of dynamical chiral symmetry 
breaking, where the pion emerged as a Goldstone bo¬ 
son and the nucleon was the fundamental degree of free¬ 
dom [28, 29]. The NJL model has subsequently been 
re-expressed with quarks as the fundamental constituents, 
making the relation with QCD evident. Importantly, the 
NJL model [31] preserves the fundamental symmetries 
of QCD. In particular, the generation of mass through 
the dynamical breaking of chiral symmetry is beautifully 
illustrated. In contrast, quark confinement is not auto¬ 
matically incorporated in the model. However, it has 
been shown that it can be mimicked by the introduction 
of an infrared cutoff in the proper-time regularization 
scheme [32-34]. The NJL model has a history of suc¬ 
cess in the description of numerous meson [31, 35, 36] 
and baryon [35, 36] properties, including the nucleon par- 
ton distribution functions [37-41] and electromagnetic 
form factors [27]. More recently these studies have been 
extended to the computation of the axial charges for 
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Figure 1. (Colour online) The NJL gap equation in the Hartree- 
Fock approximation, where the thin line represents the elemen¬ 
tary quark propagator, SQ^{k) = $ — m + ie, and the shaded 
circle represents the 4-fermion interaction. 


Figure 2. (Colour online) Bethe-Salpeter equation for 
antiquark-quark (meson) correlations in the NJL model using 
the random phase approximation. 


strangeness conserved /3-decays in the baryon octet [42] 
and possible insights into the solutions of long time enig¬ 
mas in QCD, such as the AI = 1/2 rule in kaon de¬ 
cays [43]. It is this wealth of achievement, together with 
the recent developments in lattice QCD that encourage us 
to test whether the model gives an accurate description 
of p-meson properties. 

In the application of the NJL model to the solution of 
the form factors of the p-meson, we use a two-flavor NJL 
Lagrangian which in the qq interaction channel reads: 


C = ip{i0 
+ I Gir 


— rn) ip 

{PpipY - 


{'ip^5T'tpy 




-hGp 




( 1 ) 


where t are the Pauli matrices representing isospin and 
TO = diag [to„, mj] is the current quark mass matrix. We 
assume rUu = md = rn. The fermion coupling G^r rep¬ 
resents the strength of the scalar (qq) and pseudoscalar 
( 9759 ) interaction channels and is responsible for the dy¬ 
namical generation of the dressed quark masses through 
the breaking of chiral symmetry. The strength of the 
vector-isoscalar and vector-isovector four fermion interac¬ 
tions is given by G^j and Gp, respectively. The explicit 
breaking of C/(l) axial symmetry is often modeled by the 
inclusion of an extra six-fermion determinant interaction 
term, which describes the 77 and ry' mass splitting [31], 
however, this is not directly related to our calculation 
so we do not consider it. We regularize the NJL inter¬ 
action through the proper-time regularization scheme, 
using an infrared cutoff (A//?) to remove unphysical decay 
thresholds for hadrons into quarks [32-34]. 

The dressed quark masses are given by the solution of 
the gap equation depicted in Fig. 1, which in the proper¬ 
time scheme reads 


M = m+^MG^ dr -(2) 

giving a dressed quark propagator of the form: 


nature of the interaction channel [27], where the reduced 
t-matrices for the tt, p and oj mesons read 




(,) = ^ 


-2iG 


p(uj) 


+ 2Gp(^) nyy(g2) 
X 


gpu 2 Gp(^) Iivv{q^) ^ 


(5) 


and the bubble diagrams are defined by 


npp ( 9 ^) = 6 * y Tip [75 S{k) 75 5'(fc -k 9 )], ( 6 ) 

nyy(9') 

9“ . 

\T^D[rS{k)Y S{k + q)\. (7) 


6i 


qPqV 

92 

_dyk 

(27r) 


The meson masses are given by the poles in the reduced 
t-matrices, that is 


l + 2G,npp(92 =to 2 ) = 0 , 

( 8 ) 

l + 2Gpnyy(9 " =to 2 ) = 0 , 

(9) 

1 + 2 Gt,j Uvv ( 9 ^ = ”4^) = 0- 

( 10 ) 


Expanding the full t-matrices about these poles gives the 
homogeneous Bethe-Salpeter vertices for the tt, p and w 
mesons: 

r; = v^75r„ r/'* = v^7''r„ r'i = v^7^ 

( 11 ) 

where the meson-quark-quark couplings read [30, 31, 36] 


= -^npp(92) 




pM ="^nyy(9 ) 




( 12 ) 

(13) 


III. RHO ELECTROMAGNETIC FORM 
FACTORS 


S{k)-^ = Ip - M + ie. (3) 

The description of mesons as 99 bound states in the 
NJL model is obtained via the Bethe-Salpeter equation 
(BSE) in the random-phase approximation, as illustrated 
in Eig. 2. The solution of the BSE in each meson channel 
is given by a two-body t-matrix that depends on the 


The electromagnetic current for a p-meson is parame¬ 
terized by three form factors and takes the form: 


9 “ 9 ^ 


jr^{p',p)= (P'+P) 


2 TO^ 


- FspiQy, 


(14) 
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Figure 3. (Colour online) Feynman diagrams representing the 
electromagnetic current for a meson in our NJL model. 


where the polarization of the incoming and outgoing p- 
meson is represented by the Lorentz indices a and /?, 
respectively, and p, is the photon polarization. From these 
form factors one can define three Sachs form factors for the 
p, namely, the charge [Gc{Q^)\\ the magnetic [Gm{Q^)] 
and quadrupole [G'q(Q^)] form factors, which read 

Gc(Q') = Fi(Q2) + |^Gq(Q2), (15) 

Gm(Q") = F,{Q^), (16) 

Gq(Q2) = Fi(Q2) + (1 + p) - FaiQ^), (17) 

where ri = -r^ and all form factors are dimensionless. 

I 4 

In our NJL model the p electromagnetic current is 
depicted in Fig. 3 and expressed by 


jpT/ip'^P) = * J 


(Fk 


X Tr 


X Tr 


(2^)4 

^ S{p'Fk)K^{p',p)S{p + k)T<^/S{k) 
Fk 




(27r)^ 


T'-/S{k-p)A^^{p',p)Sik-p')T^;^ S{k) , (18) 


where the Bethe-Salpeter vertices for the p are given in 
Eq. (11), k^{p,p') is the dressed quark-photon vertex and 
the trace is over Dirac, color and isospin indices. Follow¬ 
ing the calculations in Refs. [27, 30], we consider three 
versions of the quark-photon vertex, each of increasing 
sophistication; a pointlike quark-photon vertex; a vertex 
given by the solution of the inhomogeneous Bethe-Salpeter 
equation (illustrated in Fig. 4); and finally a quark-photon 
vertex which includes the pion cloud at the quark level 
(see Fig. 5). 

The pointlike quark-photon is simply given by 




T3 

2 


T 


(19) 


where g + ^ is the quark charge operator. Projecting 
onto flavour sectors the vertex is separated into two com¬ 
ponents: 




1 + T3 


ed ■ 


1 - Ta 


T 


( 20 ) 


where e„ and are the charges of the u and d quarks, 
respectively. 



y y 

= ww/viQ + ww^ • 


Figure 4. (Colour online) Inhomogeneous BSE whose solution 
gives the quark-photon vertex, represented as the large shaded 
oval. The small circle depicts the pointlike quark-photon 
driving term of Eq. (19), whereas the shaded circles with the 
double line represents the vector meson t-matrices. 


In general the quark-photon vertex is dressed by qq 
interactions in the vector channel and in the NJL model 
this dressing is described by the corresponding inhomo¬ 
geneous Bethe-Salpeter equation (see Fig. 4). From the 
NJL Lagrangian of Eq. (1) the contributions to this ver¬ 
tex come from the neutral vector mesons (p° and ui). In 
the on-shell approximation for the external quarks, the 
solution of the inhomogeneous Bethe-Salpeter equation 
of Fig. 4 is 






where the dressed quark form factors are 
.. 1 

-Ti 


yyq) 1 + 2G, 


(p) nvy(g2) 


( 21 ) 


( 22 ) 


Note, with the Lagrangian of Eq. (1) the inhomogeneous 
Bethe-Salpeter equation does not generate a Pauli form 
factor for the dressed quarks. Again projecting into 
flavour sectors gives 






T, 

(23) 


where the dressed quark form factors are given by [27] 

F^lTiQ^) = IfUQ") + If^FQ^), (24) 

F^h^iQ^) = IfUQ^) - iFip(Q2). (25) 

Finally we include pion loop corrections to the quark- 
photon vertex, as illustrated in Fig. 5, which give a vertex 
of the form 


A>^{p,p') = A[)(p,p') -b A^(p,p') 
where the flavour sector vertices {Q = U, D) read 


(26) 


^Q^P') =yFiQ{q^) 


2M 


F2q(Q2). (27) 


Note that the pion cloud generates a Pauli form factor 
for the dressed quarks and that in obtaining Eq. (27) we 
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Figure 5. (Colour online) Pion cloud contribution to the quark- 
photon vertex. The quark-photon interaction in the first two 
diagrams, represented by the shaded ovals, is given by the 
solution of the inhomogeneous BSE. The last diagram also 
includes the pion form factor determined without the pion 
cloud on the associated dressed quarks. 


have assumed the external quark lines are on-shell. The 
dressed quark form factors now read [27] 


Fiu = 


Fiuj + 


+ [A’lo; — Fip] 

n 

+ Fipfi , 







(28) 

Fid = 

-z[l 

Fiuj — 


+ [Fluj + Fip] 

n 

- Fipf^, 







(29) 

F 2 U = 

= [Fi^ 

-Fip] 

/! + 

FipH, 


(30) 

F 2 D = 

= [Fi^ 

+ Fip] 

/I- 

Fipfh 


(31) 


where for clarity we have dropped the explicit depen¬ 
dence. The renormalization factor Z is given by 


Z = 1-b 



(32) 


where S(p) is the self-energy from the pion cloud on a 
dressed quark: 



d'^k 


75 TiS{p- k) 75 Ti T„{k). 


(33) 


Here the pion propagator is approximated by its pole 
form 




i Z-^ 

— rri^ + ie 


(34) 


The contributions of the pion cloud to the quark- 
photon vertex are contained in the functions /f (Q^) and 
(i = 1, 2) of Eqs. (28)-(31). These body form 
factors are associated with the second and third diagrams 
in Fig. 5, which are respectively expressed as 

= r /|( 9 "), (35) 

+ (36) 

The analytic expressions read 

Aq^''(p',p) “ k)Y iS{p - k)-f^T^{k), 

(37) 


M 

Am 

Auv 


Gp 


Z-K 


Zuj 

0.4 

0.24 

0.645 

19.04 

11.04 

10.41 

17.85 

6.96 

6.63 


Table I. Parameters of the model together with the effective 
couplings computed from Eqs. (12)-(13). The masses are in 
units of GeV, the Lagrangian couplings in units of GeV“^ and 
the effective couplings are dimensionless. 


aW''(p',p) = Z-1 (p' +p)'^FiPL)(^2) 

X J -k)T^{p-k)-fziS{k), (38) 

where is the pion form factor determined with 

a pointlike quark-photon vertex. 

For the full calculation of the p-meson form factors 
we use Fq. (18) and the quark-photon vertex given in 
Fq. (26). For the p'^ form factors this gives 


F,pAQ^) = [FiuiQ^) - FMQ^)] fYiQ^) 

+ [F2u{QY - F2d{QY] fl(QY, (39) 

where i = 1, 2, 3 indicates each of the three form factors 
of Fq. (14). The body form factors fY are associated 
with the vector part ( 7 ^) of the quark-photon vertex in 
Fq. (27), while //’ are the body form factors associated 
with the tensor coupling (^^- 5 ^) in Fq. (27). To ob¬ 
tain the p-meson form factors that result only from the 
inhomogeneous BSF quark-photon vertex we then sim¬ 
ply set Z = 1 and the pion cloud contributions (/f(Q^), 
etc) to zero. Finally, the p form factors for a 
pointlike quark-photon vertex are then obtained by setting 
= Fip = 1. Note, all loop integrals are regularized 
using the proper-time scheme, with both an infrared and 
ultraviolet cutoff, except those of Fqs. (33), (37) and (38), 
where we take the infrared cutoff (A//j) to zero as the 
pion should not be confined. 


IV. RESULTS 

The parameters of our model are the dressed quark 
mass M; the regularization cutoffs Ajjv and Am; and the 
Lagrangian couplings Gtt, Gp and Gc^. For consistency 
with previous work we set M = 0.4 GeV (in the physical 
limit: = 140 MeV) and Am = 0.24 GeV [27, 30, 42]. 

The ultraviolet cutoff Ajjy is fit to the physical value of the 
pion decay constant and the couplings G,r, Gp, and G^i are 
fit to the physical masses of the tt, p and w mesons using 
Fqs. (8)-(10). The values of these parameters, together 
with the quark-meson couplings of Fqs. (12)-(13), are 
given in Tab. I. 

Our purpose here is to compare results within this NJL 
model with other calculations, for example, constituent 
quark models [6-14], QGD sum rules [15-18], Dyson- 
Schwinger equations [19-22] and the recent lattice QGD 
studies [25, 26]. We first focus on static electromagnetic 
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Reference 


Pp (pn) 

Qp (fm") 

This work 

0.67 

3.14 

-0.070 

Garcia Gudino [6] 

- 

2.6(6) 

- 

Cardarelli [8] 

0.35 

2.76 

-0.024 

De Melo [9] 

0.37 

2.61 

-0.052 

Melikhov [11] 

0.33 

2.87 

-0.031 

Jaus [12] 

- 

2.23 

-0.022 

Choi [13] 

- 

2.34 

-0.028 

Biernat [14] 

- 

2.68 

-0.027 

Samsonov [16] 

- 

2.4(4) 

- 

Aliev [18] 

- 

2.8(6) 

- 

Hawes [19] 

0.37 

3.28 

-0.055 

Bhagwat [20] 

0.54 

2.54 

-0.026 

Roberts [21] 

0.31 

2.14 

-0.037 

Pitschmann [22] 

- 

2.13 

- 

Owen [25] 

0.670(68) 

2.613(97) 

-0.0452(61) 

Shultz [26] 

0.30(6) 

2.00(9) 

-0.020(4) 


Table II. Comparison of the charge radius, magnetic mo¬ 
ment and quadrupole moment for various theoretical ap¬ 
proaches: phenomenological models [6], constituent quark 
models [8, 9, 11-14], QCD sum rules [16, 18], Dyson Schwinger 
equations [19-22] and lattice QCD [25, 26]. The lightest pion 
mass used in the lattice calculation in Ref. [25] is = 0.026 
GeV^, whereas for Ref. [26] it is = 0.49 GeV^. 


quantities for the meson and consider the magnetic 
moment (p.p), quadrupole moment (Qp) and rms charge 
radius ((r^)). These observables are defined by the Sachs 
form factors given in Eqs. (15)-(17), where the magnetic 
moment in nuclear magnetons (/i^v) is given by Pp = 
with Mjv the physical nucleon mass and (for 
comparison with lattice data) rup is the p mass evaluated 
at a particular pion mass; the quadrupole moment in 
units of e/Wp is given by Qp = G'q(O); and finally the 
charge radius is defined by 


(-S) = -6 


dGcm 

dQ^ 


Q2=0 


(40) 


In Tab. II we summarize results for the (r^), Pp and Qp 
of the p'^ from various theoretical approaches, together 
with our calculations using the most sophisticated quark- 
photon vertex of Eq. (26) (BSE -|- pion cloud). In general 
including the dressing of the quark-photon vertex by the 
BSE and the pion cloud increases the magnitude of /ip by 
24%, Qp by 22% and (r^) by 16% [27]. 

In comparing our results with lattice QCD we focus 
on the lattice simulation from Ref. [25], as they extend 
to the lightest pion mass, namely, uItt = 161 MeV. Our 
computations as functions of are performed by keep¬ 
ing the regularization parameters {Ajji and Ajjv) and 
the couplings {G-^, Gp and Gui) fixed, and varying the 
current quark mass that enters the gap equation. Results 
for the p mass as a function of (or equivalently the 
current quark mass) are presented in Fig. 6, where we 



Figure 6. (Colour online) NJL model results for the p-meson 
mass versus m^. Comparison is made with lattice results from 
Ref. [25]. 


find remarkable agreement between our NJL calculation 
and the lattice results of Ref. [25] . 

At the physical pion mass our values for (r^) (see 
Tab. II) differ significantly from the constituent quark 
models, one of the Dyson-Schwinger calculations and the 
result quoted in the lattice QCD computation of Ref. [26] . 
Better agreement is seen with the Dyson-Schwinger equa¬ 
tion calculation of Ref. [20]. Our result for (r^) is however 
very similar to the lattice QCD value obtained in Ref. [25] 
for a pion mass of around 161 MeV. We see that in Fig. 7 
their (r^) lies around 0.67 fm^, possibly reaching 0.7 fm^ 
in the physical limit. On the other hand the lattice QCD 
simulation of Ref. [26] uses a very large pion mass of 700 



Figure 7. (Colour online) The squared charge radius (r^) for 
the meson computed using the three levels of sophistica¬ 
tion for the quark-photon vertex: pointlike (PL), using the 
inhomogeneous BSE (BSE) and including the pion cloud (BSE 
-I- pion cloud). Comparison is made with lattice results from 
Ref. [25]. 
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Figure 8. (Colour online) The p"*" magnetic moment (pp) 
versus The curves have the same conventions as Fig. 7. 



Figure 9. (Colour online) The p^ quadrupole moment (Qp) 
versus The curves have the same conventions as Fig. 7. 


MeV, which explains its lower value for (r^), evident from 
the dependence of the lattice points in Fig. 7. The 
dependence of (r^) on in our NJL calculation, once 
the inhomogeneous BSE and pion cloud contributions 
have been included, shows remarkable agreement with the 
lattice results of Ref. [25] . One sees that the pion cloud 
contributions have become negligible for > 0.4GeV^. 

For the p"*" magnetic moment (pp) the values obtained 
by the constituent quark models are consistently smaller 
than our result of Pp = 3.14 pw, the closest being pp = 
2.87fiN from Ref. [11]. The earlier Dyson-Schwinger 
equation study in Ref. [19] shows good agreement with 
our work. For the lattice simulation of Ref. [25] the 
discrepancy with our result is sizeable near the physical 
limit. However, the evolution of our result with shown 
in Fig. 8 is in good agreement with the lattice QCD 
calculations except at their lightest pion mass. Again, as 
in the case of (r^), the effect of the large m,r in Ref. [26] 





ml (GeV2) 

Figure 10. (Colour online) The p^ Sachs form factors as a 
functions of ml at — 0.16 GeV^. The curves have the 
same conventions as Fig. 7. 


is to produce a small value of pp, as evident from Fig. 8. 

Finally we find a large quadrupole moment comparable 
to the Dyson-Schwinger equation results of Roberts et 
al. [21] and Hawes et al. [19]. The lattice QGD result of 
Ref. [25] is approximately ~30% smaller than our result 
near the physics point, as illustrated in Fig. 9. However, 
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(GeV2) 

Figure 11. (Colour online) The Sachs form factors for 
= 0.49 GeV^. Comparison is made with lattice results 
from Ref. [26] and the curves have the same conventions as 
Fig. 7. 


as a hypothesis for the difference we suggest that it may be 
worthwhile to investigate the effect of the lack of spherical 
symmetry on the lattice simulation, considering that the 
quadrupole moment reflects the shape of the p. 

Comparison with the lattice simulation of Ref. [25] for 
the evolution of the p'^ Sachs form factors with mj, at 
a fixed = 0.16 GeV^, is made in Figs. 10. The charge 
form factor, Gc, is in good agreement with the lattice 
QCD points, when both the inhomogeneous BSE and 
pion cloud dressing are included. On the other hand, for 
the magnetic form factor Gm, the BSE results alone have 
better agreement with lattice and the pion cloud causes 
an overestimate. The deviations are still small however, 
considering the simplicity of the calculation. The devia¬ 
tion from the lattice simulation data for Gq is possibly 
explained by the same reason behind the disagreement 
with Qp, that is, the lack of spherical symmetry in the 
lattice simulation. 

A final comparison is made in Fig. 11 for the Sachs form 
factors as a function of for a pion mass of — 0.49 
GeV^, where the lattice results are from Ref. [26]. We 
find that our model qualitatively describes the p^ form 
factors obtained from the lattice computation. Once again 
the addition of the pion cloud causes an overestimate of 
Gm{Q^) and the magnitude of Gq((5 ^) also appears too 


large. 

V. CONCLUSIONS 

We computed the electromagnetic form factors of the 
meson using an NJL model that simulates aspects of 
quark confinement. The quark-photon vertex is studied in 
three levels of sophistication: pointlike dressed quark, via 
the inhomogeneous BSE and also including corrections 
from a pion cloud. The results are qualitatively in good 
agreement with the recent lattice QCD computations. 

The main level of disagreement comes from the 
quadrupole moment and the corresponding form factor. 
We suggest that lattice QCD studies of this type should 
look at the possible effects of the lack of spherical sym¬ 
metry of a cubic lattice in the quadrupole moments and 
form factors. It would certainly be helpful to have further 
lattice studies over a range of pion masses and momen¬ 
tum transfers. Experimental measurements would also 
be extremely valuable. 

Therefore, the present work on the p-meson structure 
and the progress in the computation of the electromag¬ 
netic form factors of the tt and K, including the pion 
cloud, reported in Ref. [.30], support the importance of 
the model as a tool to describe hadronic structure. In 
addition, the NJL model is a quantum field theory where 
calculations are relatively straightforward and it gives 
good results when compared to more sophisticated meth¬ 
ods that require much more resources, such as lattice 
QCD. These advantages are useful in order to perform 
larger calculations in problems such as the description of 
hadrons in the nuclear medium, as required, for example, 
to explore the properties of neutron stars. In such cases 
the NJL model serves as a very useful tool to guide possi¬ 
ble future computations of lattice QCD and other more 
sophisticated approaches. 
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